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Discrete energy levels of ultrasmall metallic grains are extracted in single-electron-tunneling- 
spectroscopy experiments. We study the response of these energy levels to an external magnetic 
field in the presence of both spin-orbit scattering and pairing correlations. In particular, we investi¬ 
gate g-factors and level curvatures that parametrize, respectively, the linear and quadratic terms in 
the magnetic-field dependence of the many-particle energy levels of the grain. Both of these quan¬ 
tities exhibit level-to-level fluctuations in the presence of spin-orbit scattering. We show that the 
distribution of g-factors is not affected by the pairing interaction and that the distribution of level 
curvatures is sensitive to pairing correlations even in the smallest grains in which the pairing gap is 
smaller than the mean single-particle level spacing. We propose the level curvature in a magnetic 
field as a tool to probe pairing correlations in tunneling spectroscopy experiments. 
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I. INTRODUCTION 

Discrete electronic energy levels in individual nano¬ 
scale metallic grains (nanoparticles) can be measured by 
means of single-electron-tunneling spectroscopyji Since 
the pioneering experiments in aluminum grains spec¬ 
tra of nanoparticles made of a variety of materials were 
measured. In particular, the magnetic-field response of 
energy levels of grains with spin-orbit scattering has been 
studied in a number of experiments 

Spin-orbit scattering breaks spin symmetry. However, 
time-reversal symmetry, leading to Kramers degeneracy 
of the single-electron levels in the absence of a magnetic 
field, remains a good symmetry. Since atomic-scale ir¬ 
regularities of a typical grain destroy all possible orbital 
symmetries, there are generally no additional degenera¬ 
cies, and the single-electron levels form doublets. An ex¬ 
ternal magnetic field breaks time-reversal symmetry and 
leads to splitting of the doublets. The energies of the 
upward- (efe+) and downward-moving (cfc-) levels of a 
doublet Efe in a weak magnetic field B are parametrized 
by the g-factor gk and level curvature Kk (at zero field) 
as 

£k±{B) = efc(O) ± -QkfJ-BB -\- -HkB'^ A 0{B^), (1) 

where /is is the Bohr magneton. 

When spin-orbit coupling is negligible, the spin is a 
good quantum number, orbital magnetism does not con¬ 
tribute to g-factorS )^°d^ and = 2 for all levels. How¬ 
ever, in the presence of spin-orbit coupling, the g-factor 
exhibits level-to-level fluctuations^ip— and depends on 
the magnetic-field directionJ^d^ In addition, the spin 
contribution is suppressed and may become comparable 
to the contribution from orbital magnetism,— which 
is finite when spin symmetry is broken. The statistics of 


single-particle g-factors have been studied using random- 
matrix theory (RMT)fii“— and are generally in good 
agreement with spectroscopy experiments on noble-metal 
nanoparticles»i“— However, the understanding of the av¬ 
erage value of the g-factor is still lacking. 

Another signature of spin-orbit scattering is the non¬ 
linear second-order correction (in the magnetic field) to 
energies, which results in avoided crossing of energy lev- 
elsi^d^ Although this correction is not identically zero 
in the absence of spin-orbit scattering because of orbital 
magnetism, it becomes appreciable only when spin-orbit 
scattering breaks spin-rotation symmetry, leading to fi¬ 
nite matrix elements of the spin operator between states 
of different doublets. The distribution of level curvatures 
was measured in gold nanoparticles^ and found to be in 
agreement with RMT predictions— “— 

The above single-particle level results are valid for 
grains that are well described by the constant-interaction 
(Cl) model, in which the electron-electron interaction 
is described by the classical charging energy e^A^/2C', 
where N is the electron number and C is the capaci¬ 
tance of the grain. In the Cl model, the energies ex¬ 
tracted in a tunneling-spectroscopy experiment reduce to 
the single-electron energies, as is the case of noble-metal 
nanoparticles— ^ However, when the interaction effects 
beyond the Cl model are important (e.g., in supercon¬ 
ducting or ferromagnetic materials), the extracted en¬ 
ergies in the tunneling-spectroscopy experiments do not 
reduce to single-particle quantities. 

Interaction effects in a chaotic or weakly disordered 
grain are described by the universal Hamiltonia n^^d^ in 
the limit of large dimensionless Thouless conductance 
gxh = Etii/S 1. Here Axh is the Thouless energy, 
determined by the time it takes for an electron to cross 
the grain, and S is the mean spacing between single¬ 
particle Kramers doublets. The one-body part of the 
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universal Hamiltonian follows RMT statistics;^ while 
the interaction consists of universal terms that are con¬ 
sistent with the symmetries of the one-body Hamilto- 
niani ^^d^i^^ In the absence of spin-orbit scattering and 
orbital magnetic field, these terms are the charging en¬ 
ergy (as in the Cl model), pairing interaction, and fer¬ 
romagnetic spin-exchange interaction. The mesoscopic 
transport properties of a grain described by this univer¬ 
sal Hamiltonian were studied in Ref. [12, and its thermo¬ 
dynamic observables were calculated in Ref. [12 using a 
Hubbard-Stratononvicb24i2^ decomposition for the pair¬ 
ing interaction and a spin projection method2£ for the ex¬ 
change interaction. In the absence of pairing correlations 
(i..e, when only charging energy and spin-exchange cor¬ 
relations are present), various observables can be calcu¬ 
lated in closed form by using the Hubbard-Stratonovich 
transformation employed in Refs. [12 and [12 or by the 
spin projection method. 

In the presence of strong spin-orbit scattering, spin- 
rotation symmetry is completely broken, and the ex¬ 
change term is absent from the universal Hamiltonian. 
The effect of the exchange interaction on the ^-factor 
statistics in the crossover between weak and strong spin- 
orbit scattering was studied in Refs. [12 and Large 
^-factors were recently observed in ferromagnetic cobalt 
nanoparticles 

Spin-orbit scattering preserves time-reversal symme¬ 
try and therefore does not suppress pairing correlations. 
In this work we explore the effect of an attractive pair¬ 
ing interaction on the magnetic-field response of many- 
particle energy levels of a grain in the presence of strong 
spin-orbit scattering, where the exchange interaction is 
negligible. 

We assume the one-bottleneck geometr y^i^^i^° of the 
tunneling-spectroscopy experiments, for which the rate 
of tunneling onto the grain is much smaller than the 
rate of tunneling off the grain. In this limit, the cur¬ 
rent is determined by processes of tunneling onto the 
grain in its ground state. When the tunneling occurs 
onto a grain with even particle number N^, the differ¬ 
ence AEq^o = Eq‘~^^ — between the energies of a 
twofold-degenerate level H for an odd electron number 
A^e + 1 and the non-degenerate even ground state E^'^ 
is measured. This energy difference splits in an exter¬ 
nal magnetic field. We define the many-body ^-factor g 
and level curvature k by generalizing the single-particle 
expression in m 

AEa,o{B) = ASo,o(0) ± + 0{B^). (2) 

In the absence of pairing correlations, the TVe-particle 
ground state consists of doubly occupied orbitals up to 
the Fermi level, while the Ne + 1-particle level has one 
extra singly occupied orbital fco above the Fermi energy 
as is illustrated in panels (a) and (c) of Fig. [T] In this 
limit, the (^-factor and level curvature defined in Eq. © 
reduce to the single-particle values gk^ and Kkg of the 
singly occupied orbital fep. In the presence of pairing 


correlations, the single-particle occupation numbers are 
no longer good quantum numbers. The even ground state 
becomes a superposition of Slater determinants describ¬ 
ing fully paired noninteracting states as shown schemat¬ 
ically in Fig.[T])b). A final odd-particle-number state is a 
superposition of states of good occupation numbers with 
one singly occupied “blocked” orbital /cq, which is the 
same in all the Slater determinants of the superposition^^ 
as shown in Fig. [T^d). 

Here we will show that the (/-factor as defined in ([5]) is 
identical to the single-particle (/-factor (/^^ of the blocked 
orbital /cq • On the other hand, the level curvature differs 
from its single-particle value, and its statistics are highly 
sensitive to pairing correlations. The main origin of this 
sensitivity is the change of the density of states induced 
by pairing correlations. 

Another motivation of our studies is to identify sig¬ 
natures of pairing correlations in ultrasmall supercon¬ 
ducting grains whose single-particle mean level spacing 
5 is comparable to or larger than the bulk pairing gap 
A of the grain. Anderson’s criterion^^ states that su¬ 
perconductivity is no longer possible once A/5 < 1, 
and the tunneling-spectroscopy experiments showed no 
traces of an excitation gap in the smallest aluminum 
grains that satisfy this conditioni^ This criterion signifies 
the breakdown of the Bardeen-Cooper-Schrieffer (BCS) 
mean-field theory^l in the fluctuation-dominated regime 
A/d < I. A better probe to detect pairing correlations 
in this regime are thermodynamic observables such as 
the heat capacity and spin susceptibility j^^i^^'^^ but they 
are difficult to measure experimentally. Here we pro¬ 
pose the level curvature as a sensitive observable to probe 
pairing correlations in grains with spin-orbit scattering. 
Its advantage is that it can be measured directly in the 
tunneling-spectroscopy experiments. 

The outline of the paper is as follows. In Sec. m we 
introduce our model and discuss the many-body eigen¬ 
states of its corresponding Hamiltonian. In Sec. m we 
prove the robustness of (/-factor statistics with respect 
to pairing correlations. In Sec. IIVI we calculate the 
level curvature and discuss its statistics. We conclude 
in Sec. [Vj 

II. MODEL 
A. Hamiltonian 

The universal Hamiltonian in the presence of strong 
spin-orbit scattering has the following form for fixed par¬ 
ticle number 

H = Y,ek (c^iCfci + cl^ck2) - GP^P , (3) 

k 

where 

= ^ A 1 A 2 and P = Y1 Ck2Cki (4) 

k k 
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G = 0 

(a) _ (b) 


G/O 


Cl + C2 + C3 + ... 


(c) _ (d) _ _ 

ko ko [- 5 ——^ ^ ^ ~ ^ 

c; +c^^*^+c^ 


FIG. 1. Top: schematic description of the ground state of 
the Hamiltonian for an even number of particles in (a) 
the absence of pairing (G = 0) and (b) the presence of pair¬ 
ing (G 7 ^ 0). Bottom: a possible odd-particle-number state 
after the tunneling of an additional electron onto the even 
ground state in (c) the absence and (d) the presence of pair¬ 
ing. The solid lines describe doubly degenerate single-particle 
levels and the solid circles denote particles occupying these 
levels. 


are pair creation and annihilation operators. The single¬ 
particle states appear in degenerate time-reversed pairs 
|fcl)i and \k2)i with single-particle energy e^. When spin 
symmetry is broken, they are no longer eigenstates of the 
spin-projection operator. We adopt the phase convention 

f |fcl)i = \k2)i and f |A:2)i = -|fcl)i , (5) 

where T is the time-reversal operator. The interaction 
term in Eq. describes the scattering of electron pairs 
between such Kramers-degenerate orbitals and leads to 
BCS superconductivity in the bulk limit A/5 ^ 1. The 
exchange interaction is suppressed for strong spin-orbit 
scattering and is absent in Eq. ©■ 

The coupling of electrons in the grain to an external 
magnetic field along the z-axis is described by the Zee- 
man term 


V{B) = -M,B , (6) 

where Mz is the z-component of the magnetic-moment 
operator of those electrons. In general, the response of 
levels in the grain to an external magnetic field depends 
on its direction^i^iH but we do not study this effect here. 

In a chaotic or a weakly disordered grain, the single¬ 
particle levels follow RMT statistics. In the absence 
of an external magnetic field, time-reversal symmetry is 
preserved, and the single-particle spectrum is described 
by an ensemble that interpolates between the Gaussian 
orthogonal ensemble (GOE) and the Gaussian symplec- 
tic ensemble (GSE) In the absence of spin-orbit 

scattering, spin-rotation symmetry is preserved and the 
corresponding ensemble is the GOE. In the presence of 
strong spin-orbit scattering, spin-rotation symmetry is 
completely broken and the corresponding ensemble is the 


GSE. We note, however, that most of our qualitative con¬ 
clusions are independent of the particular statistics of the 
single-particle eigenstates and eigenenergies. 

In our quantitative calculations of level curvatures in 
Sec.EVl we make the following additional assumptions: 
(i) spin-orbit scattering is strong, i.e., the single-particle 
Hamiltonian follows the GSE statistics, (ii) the orbital 
contribution to the magnetic moment is negligible, i.e., 
Mz = 2 ij,bSz, where Sz is the z-component of the spin 
operator, and (iii) in a finite-size model space with Nsp 
Kramers-degenerate orbitals, the pairing coupling con¬ 
stant is determined from the relatio n^^i^^ 


arcsinh ^ ^ 


B. Many-body eigenstates 

In the absence of pairing correlations, each many- 
particle eigenstate of the Hamiltonian m is a state 
with well-defined occupation numbers of single-particle 
orbitals. We will denote such noninteracting state by 
\m, kiai ... ktat)N, where m = {mi,..., m^} is a set of r 
doubly occupied orbitals, ki,... ,kt are t singly occupied 
orbitals {ki ^ m) with corresponding labels ai,..., at 
{at = 1,2) distinguishing Kramers-degenerate single¬ 
particle states, and N = 2r + t is the total number of 
particles. 

In the presence of pairing correlations, each eigenstate 
is a superposition of these noninteracting states. Since 
pair scattering cannot affect the singly occupied orbitals, 
the orbitals fci, ai,..., fct, a* are good quantum numbers 
and are the same in all the noninteracting states com¬ 
prising the superposition of any given eigenstate of the 
Hamiltonian ([3]) . This is a manifestation of the blocking 
effect of the pairing interaction.^ An A-particle eigen¬ 
state can thus be written as 

, kiai ... hat) N = ^ C%^\m,kiai.. .ktat)N ■ 

TTLiki^rh 

( 8 ) 

The sum in Eq. (|S]) runs over all sets m of {N — t)j2 
doubly occupied orbitals that do not contain any of the 
singly occupied orbitals h. The coefficients C? can be 
chosen to be real. The energy of the state (|5]) and the cor¬ 
responding coefficients C? are independent of ai. Since 
each ai can assume two values, the corresponding many- 
particle level is 2*-fold degenerate. 

The ground state |0)jy of © for an even number of 
particles is described by a superposition of fully paired 
noninteracting states with no singly occupied orbitals 

|0)^. (9) 

m 

and is schematically shown in Fig. [IJb). A possible odd 
state after tunneling onto the grain in |0)^^ 
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must have non-zero overlap with |0)^ for at least 
one single-particle orbital \ka)i. This, in combination 
with the blocking effect, dictates that such an eigenstate 
has exactly one blocked (i.e., singly occupied) orbital as 
illustrated in Fig. [T])d). These states form doublets and 
can be expressed as 

I'k, fcoao)Afe-i-i = ^ C^\rh,koao)N^+i ■ (10) 

rn-.ko^rh 


III. p-FACTOR 


The many-particle g-factor defined in Eq. ([2]) depends, 
in general, on the linear corrections to both the odd- 
and even-particle energies. However, since the even state 
|0)Ar,, is invariant under time reversal and is odd. 





( 11 ) 


and the linear correction to the even ground-state en¬ 
ergy Eq‘ vanishes. The g-factor that corresponds to the 
|0)^^ —>■ |fl)jv,,+i = transition is then deter¬ 

mined by the 2x2 matrix of Mz written for the doublet 
jdt, From Eq. (flUl) . we find 


koa ^Mz 


dt, koa' 


JVe-l-l 


^iTi^rh 


m, koa 


Mz 


fh', koa' 


Ne + l 

( 12 ) 


Using 

-^2 = ^kp,k>p>^\p^k'p' ) (13) 

kp,k' p' 

where are the single-particle matrix elements of 

Mz: we have 


TO, koa 


Mz m',koa'} =5^ 
Ne + l 


M, 


koocykoOi' 


According to Eq. dm), the matrix elements of the 
magnetic-moment operator between any two many- 
particle eigenstates that belong to the same Kramers 
doublet reduce to the single-particle matrix elements 
of the blocked orbital. Therefore, the g-factor for the 
|0)Ar^ —>■ |4>, fcoa)Ar^+i transition is exactly the single¬ 
particle g-factor of the blocked orbital feg. 

The above result is quite general. It is independent 
of the statistics of the single-particle levels, strength of 
spin-orbit scattering, or relative weights of the spin and 
orbital parts in the magnetic moment. It follows from 
the blocking effect of the pairing interaction and from 
time-reversal symmetry. The blocking effect makes one 
orbital ko special in the odd state and separates the mag¬ 
netic moment of that state into the contributions from a 
single electron on the orbital ko and from the remaining 
fully-paired electrons. Time-reversal symmetry makes 
the contributions from the paired electrons to the mag¬ 
netic moments of the odd and even states zero. 

Our result can be verified in a tunneling-spectroscopy 
experiment by measuring (^-factors at different particle 
numbers (achieved by varying the gate voltage). Since 
the g-factor is completely determined by the blocked or¬ 
bital ko of the final state of a tunneling process |0)Ar^ —>■ 
jdt, fcoa)Ar^+i, it remains the same for any other allowed 
process |0 )jy^ —>■ |^,fcoa)jv„+i with Ne ^ A^e as long as 
the blocked orbital of the final state is the same. 

Another consequence of our result is that the measured 
distribution of the g-factor reduces to the distribution of 
the single-particle g-factor. For sufficiently strong spin- 
orbit scattering, the latter distribution has a universal 
formiiii^ that should not be affected by pairing correla¬ 
tions. We note that the RMT predicts the distribution 
of g-factors measured in units of their average values (g ), 
but not the average values themselves. 

A more general pairing interaction with 
orbital-dependent coupling constants of the form 
Sfc.fc' Gkk'cl^cl^Ck' 2 Ck'i leads to blocking effects as 
well and does not modify the g-factor. Therefore, 
the g-factor can be used to probe electron-electron 
correlations beyond any such generalized pairing model. 


+ Saa' ^ 

m^rh 


(14) 


Since the single-particle states |to 1) and |to 2) are mapped 
onto each other under time reversal [see Eq. (O] and Mz 
is odd under time reversal. 


,ml ,m2 


= 0 . 


(15) 


We conclude 


4', koa 


Mz 


dt, koa' 


We-l-l 


= M, 


kQOi,kQOL' 5 


(16) 


where we have used the normalization condition 


In the absence of spin-orbit scattering, the exchange 

/K 2 --V 

interaction —JsS (where S is the total spin of the uni¬ 
versal Hamiltonian) is consistent with the symmetries of 
the single-particle Hamiltonian and does not destroy the 
blocking effect. Finite spin-orbit scattering makes the 
structure of this term complex when written in the basis 
diagonalizing single-particle Hamiltonian [i.e., the |fca)i 
basis in Eq. (I3])[ and destroys the blocking effect. The 
many-body eigenstates become superpositions of nonin¬ 
teracting states with different singly occupied orbitals, 
so the magnetic moment of such a state can no longer 
be reduced to a single-particle magnetic moment. The 
g-factors and their distribution are affected by the ex¬ 
change interaction, as was studied in detail in Refs. 
and Isd. 
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IV. LEVEL CURVATURE 


Using second-order perturbation theory, the level cur¬ 
vature for the transition | 0 ) 7 Ve |U)Ar^+i is given by 


K = K 


N^ + l 
n 


^Ne 


- Kr, 


(17) 


where 


and 



/ 


\ 

^JVe + l _ OV"'. 

(u 


U') 

/We + l 


^ / ■' p-Ve + l p-^e-t-1 


N, 


= 2E' 


(o 

Mz 

0') 

\ 


/ Ne 


j^Ne 

-C/Q ^0' 


(18) 


(19) 


Here the sums run over the many-electron eigenstates 
with energies different from or Eq". The expres¬ 

sions m and m can be thought of as the curvatures 
of the odd and even states, respectively. 

In the G —>■ 0 limit, the curvature (HZl) reduces to the 
single-particle curvature 


Kko 


k^ko 


\M, 


kol^kl 




Cfco ~ 


( 20 ) 


where fcp is the blocked orbital in the final state |U) 7 Ve-i-i- 
The corresponding distribution P^Kko) is symmetric 
around zero and was calculated analytically in the GSE 
limit i^^d^ The right (left) tail of P{Kko) is determined 
by the probability that the energy of the orbital below 
(above) /cq is very close to efcj,, which in the GSE limit re¬ 
sults in P{Kko) ^ 1/^fco ™ distribution)!^ 

In the presence of pairing, P{k) changes and depends 
on the final level |U)Ar^_|_i. When this level is the odd 
ground state | 0 ) 7 Ve-i-i, both contributions and (ED 
are negative. Therefore, n for the transition |0)Ar^ —>■ 
|0)Ar^+i is positive when and vice versa. 

Since pairing induces a gap in the even excitation spec¬ 
trum only, is suppressed by pairing, while 

is not, making the distribution asymmetric with a neg¬ 
ative median value. We verify this qualitative reasoning 
by calculating level curvatures exactly and in a gener¬ 
alized BGS approach for transitions to both the ground 
and excited states. 


A. Formalism 


(a) 




kp- 


k < kp 


(b) 


^) n^+i |n')Ar,+i 


^k'> kp ko[- 
kp - 


k < kp 


(c) 


|U)Ar,, + l |n')v. 


+ 1 


(d) 


|U);V, + 1 \^')n^ 


+ 1 


-k > kp, ^ ko 


ko ■ 
kp ■ 


ko ■ 
kp ■ 


k < kp 


FIG. 2. Pairs of many-electron states contributing to the 
second-order perturbative expressions m and m of the 
level curvature in the G = 0 limit. For an even grain, the two 
states in a pair must differ by a particle-hole excitation that 
reduces by one the number of doubly occupied levels [panel 
(a)]. For an odd grain, these states may relate to each other 
as in (a) [see panel (b)] or they can have the same number of 
doubly occupied levels but different blocked orbitals [panels 
(c) and (d)[. 


1. Noninteracting limit 


Let kp he the lowest empty orbital in the noninteract¬ 
ing even ground state. Then, the states |0')Afe contribut¬ 
ing to Eq. (|T^ are particle-hole excitations |0')Afe = 
^k'a''^ka\h)Ne with k' > kp and k < kp. A possible 
pair of states |0)Ar^ and |0')Afe is shown schematically 
in Fig. Ha). The corresponding matrix element between 
these states is (0'|Afj;|O)Ar^ = M^,^, and their energy 
difference is E^” — E^f = Ck — ek> ■ From time-reversal 
symmetry \M^,i^k 2 \ = |d^f' 2 .fcil and \M^, 2 ^k 2 \ = |d^fe'i.fcil 
so the curvature of the even state is 


K 


Ne 

0 


G=0 


^ Mkk' 


( 21 ) 


where 


Mkk' — -I- ■ (22) 

There is a similar contribution to the curvature dT51) of 
an odd state |U) 7 Ve-i-i with a blocked orbital ko > kp [see 

Fig.Hb)] 


K 


( 1 ) 

n 


G=0 


= 4 E 

kKkp <k' ,k'^ko 


Mkk' 
Cfe — Cfe' 


(23) 


Below we demonstrate the many-body formalism for 
the simple noninteracting case G = 0, and then discuss 
the interacting case of finite G. 


There is also a second contribution arising from pairs of 
states that have different blocked orbitals but the same 
number of doubly occupied levels [shown in Figs. H^) 
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andlUd)]: 


3. BCS formalism 


K 


( 2 ) 

n 


G=0 


2 E 

k'^kp ,k^ko 


Mkpk 
^ko ^k 


+ 2 E 

kcikp 


Mkpk 

Cfc ~ Cfco 


(24) 


Taking the difference between the odd and even curva¬ 
tures [see Eq. (EZD], many of the terms cancel and we 
obtain Eq. ((^ . 


2. Exact formalism 


Assuming the many-electron spectrum of the pairing 
model ([3|) is known, the matrix elements of are evalu¬ 
ated in Appendix The states contributing to the even 
curvature IT^ are all the states ,kak'a')Ne with two 
blocked orbitals. The sum in the odd curvature m 
is over all the states \'i>',koaokak'a')Ne+i with three 
blocked orbitals, one of which is ko (the blocked orbital 
in |n)Ar,,+i = |4', fcoao)Ar,,+i), and over all the states 
|4>", fca)Ar^+i with a single blocked orbital k ^ ko. The 
even curvature is then given by 




= 4 EE 

k<k' ^ 


Mkk' 


Bkk’ [0, 




(25) 


'I'fefe' 


and the expressions (1^ and (IMl) for the odd curvature 
change to 


(1) . Mkk'\Bkk'['i’ 

=4 2 ^ E jrN.+i _ jrN.+i 

k<k'-.k,k'^ko 3" ^3''/cofcfc' 


(26) 


To calculate the level curvature in larger model spaces 
(up to Nsp ^ 200), we employ a generalized BCS ap- 
proach;22ii2iiL in which the blocking effect in the odd 
states is partly accounted for. The computational time 
required to solve each BCS equation and the number of 
relevant many-electron states in the sums (1181) and m 
scale as powers of Ngp, enabling calculations in much 
larger model spaces. Below we summarize this method, 
and more details are presented in Appendix |5] 

The BCS ground state for even particle number Nf, is 
given by 


|BCSe) = + ^fc4i42)|vac), 

k 

where jvac) is the vacuum state. Here, 



where 


(28) 


(29) 


Cfe = Cfe - Me - Gvl , (30) 

and Ek are the quasiparticle energies 

Ek = + A? . (31) 

The pairing gap Ag and the chemical potential /ie are 
determined from the self-consistent BCS equations 

= and Ne = (32) 

k k 


and 


^( 2 ) _ 9 Mkko\Dkko [^, 

^ 2 -^ 2 -^ ZT-ATe + l _ zT-ATg + l 


(27) 


Here, the denominators contain the differences between 
the corresponding many-body eigenenergies. The many- 
particle contributions to the matrix element ilfcfc/[0,4('] 
and Dkko [^j are given, respectively, in Eqs. (IA5|) and 
(lAlOl) of Appendix lAl When G = 0, they are identically 
zero for most of the terms, and are equal to 1 for the 
terms shown in Fig. [51 

In the noninteracting limit, the number of many- 
electron states |H')jVe-i-i and |0')Afe contributing to 
Eqs. ([T^ and (fT^ scales as a power of Asp, the num¬ 
ber of single-particle orbitals in the model space. In the 
presence of pairing, however, the number of such states 
scales combinatorially with Ngp] any state with an al¬ 
lowed configuration of the blocked orbitals has generally 
a nonzero contribution. Furthermore, the computational 
effort required to evaluate each separate many-electron 
matrix element of has a combinatorial dependence 
on Nsp as well. We have performed exact calculations 
for model spaces with Ngp <13. 


The excited states |0')Afe contributing to Eq. (fT511 are 
the two-quasiparticles excitations with the excitation en¬ 
ergies Eq‘ — Eq" given by the sums of the quasiparticle 
energies m- The curvature of the even ground state 
becomes (see Appendix IbI) 


Ne 


BCS 


- 4 E 

k<k' 


Mkk'jukVk' - Uk'VkY 
Ek + Ek' 


(33) 


where Mkk' is defined in Eq. (l22|) . Here, the many- 
electron contribution (ukVk' — Uk'VkY results in the sup¬ 
pression of the terms with far above or Ck' far below 
the Fermi level, similarly to the noninteracting restric¬ 
tion k < kp < k'. This restriction is now lifted within 
the pair-scattering energy window of ^ 2Ae around the 
Fermi level. When there is no gapped solution to the 
BCS equation (i.e., when Ag = 0), Eq. (l33)) reduces to 
the noninteracting result (I21|) . 

The result (1331) is very similar to Belyaev formula^ for 
the nuclear moment of inertia and to the expression for 
the zero-temperature spin susceptibility of a supercon¬ 
ductor with spin-orbit scattering derived by Anderson4^ 
These formulas describe the suppression of the respective 
quantities caused by pairing correlations. 
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For odd particle number, we perform blocked BCS cal¬ 
culations to account for the reduction of pairing correla¬ 
tions in the odd-particle number states. We consider only 
those final states |fl) n^+i that have the lowest energy for 
a given blocked orbital fcp; the peak heights in a tunneling 
spectroscopy experiment that correspond to transitions 
to other states are suppressed in the BCS limit A/J ^ 1 
(see Appendix ICl) . 

The variational lowest-energy state with a blocked or¬ 
bital fco is given by 

|BCS, fcoao) = 4oa^ {uk^k + |vac). 

k^ko 

(34) 

The corresponding pairing gap and the chemical po¬ 
tential /Tfco are now determined from 

^ UkokVkok and Ag = 2 ^ . (35) 

k^ko k^ko 


The parameters Ukgk, Vk^k, Cfeofei and Ek^k are defined 
similarly to the even case for k ^ [Eqs. ([H]), (150]) . and 
(ISTl) above]. We approximate any other relevant state by 
a quasiparticle excitation on top of one of the states (IMt . 

Similarly to Eq. (1551) . the contribution (1^ to the odd 
curvature from pairs of states with different numbers of 
blocked orbitals becomes 


K 


( 1 ) 

n 


BCS 


-4 E 


^kk' {‘^kok'^kok* '^kok''^kok) 
Ekgk + Ekok' 


(36) 

The second contribution (EH) reduces to (see Ap¬ 
pendix |B]) 


.(2) 


BCS 


-2 E ^^0 

k^ko 


{'^kok'^kko '^kok'^kko') 


('^kok'^kko Ukko‘k^kok') 
ABCSABCS 


(37) 


where 


Ek^^ — E (^^ofc “ Ekok + Gvf.^k) + + k-koEle + ^ko 

k^ko 

(38) 

is the BCS energy of the state (l34l) . For each blocked 
orbital k ko, two odd-particle-number doublets con¬ 
tribute to Eq. (EH). One of them is the lowest-energy 
level with the blocked orbital k, and the other is a two- 
qusiparticles excitation on top of |BCS, ka) in which the 
two quasiparticles occupy the same orbital k^, so the 
number of blocked orbitals does not change. If the block¬ 
ing effect in the BCS calculations were ignored, the de¬ 
nominators of these two contributions would be Ek — Ekg 
and Ek Eko , and the entire expression would resem¬ 
ble the odd-particle-number nuclear moment of inertia 
at zero temperature derived in Ref. H . 


B. Numerical simulations 

We carried out both exact and BCS calculations of 
the level curvature. For each calculation, we generate 
an ensemble of 1000 2Ai.i„t x GSE random matri¬ 

ces with Aj-mt = 201 degenerate eigenvalues each. The 
matrices are diagonalized by symplectic transformations 
using the phase convention ([S]) for the eigenvectors. For 
each matrix, we form the model space of N^-p < Armt 
single-particle Kramers doublets by taking Ngp levels in 
the middle of the spectrum and unfolding^ the eigen¬ 
values. We use Asp = 13 in the exact numerical cal¬ 
culations and Asp = 121 in the BCS calculations. As¬ 
suming the orbital contribution to the magnetization is 
negligible, we use the eigenvectors to calculate the single¬ 
particle matrix elements of = 2^,bSz in the basis of 
2Asp states forming the chosen doublets. We perform the 
many-body calculations for half-filling (i.e., Ag = Asp — 1 
for odd Asp). In the exact calculations, we use the Lanc- 
zos algorithm^ to calculate the relevant many-electron 
eigenfunctions and eigenenergies. 

We select the odd eigenstates based on the 

heights of the differential-conductance peaks, which are 
calculated as discussed in Appendix [C] We select only 
those states whose peak height is at least 10 % of the 
average peak height for an allowed transition in the non¬ 
interacting limit. In that limit, a final state |n) 7 Ve+i has 
the lowest energy among all the states with the same sin¬ 
gle blocked orbital fco {ko > kp)- Similarly, in the BCS 
limit Is./5 ^ 1, the peak heights are nonzero only for the 
lowest-energy states for a given fcp, although now kp may 
also lie below kp (the peak height is suppressed when 
efcp — efco ^ A). In the BCS calculations, we consider 
only such lowest-energy states, while in the exact calcu¬ 
lations we consider general states with a single blocked 
orbital. 

We sort all resolvable odd-particle-number doublets 
according to their energies. We refer to the curvature 
for the transition | 0 )Ar^ —>■ |n)Ar^+i to a state in the n- 
th doublet as to the curvature of the n-th differential- 
conductance peak. In this convention, the first peak usu¬ 
ally corresponds to the tunneling into the odd ground 
state. 

We characterize each curvature distribution by its me¬ 
dian value Kmed and midspread (i.e., the width of the 
middle part of the distribution with 50% of the total 
area) d^, which are robust measures for a distribution 
with possibly long tails. 

To eliminate the nonphysical dependence on the 
random-matrix size, we express our results in units of 
do, the midspread of the single-particle level curvature 
distribution. In the GSE limit, this do is related to the 
(/-factor statistics byil 

do ~ 0 . 68 ^( 52 ). (39) 

This relation also holds when there is an orbital contri¬ 
bution to the magnetization. 









When the single-particle levels follow RMT statistics, 
the average matrix elements of the spin operator do not 
decay with the energy separation between two single¬ 
particle levels, and the contribution from single-particle 
levels far from kp may be important. Using BCS calcu¬ 
lations in the presence of pairing correlations, we found 
that the median values converge slowly with increasing 
values of iVgp, and we were unable to determine whether 
convergence is reached even at A^sp 200. In contrast, 
the convergence of the midspreads is clearly reached at 
iVgp ^ 100, and the results at iVgp ^ 10 are already quite 
close to their asymptotic values. We, therefore, consider 
the median values only as qualitatively correct in both 
the exact and the BCS calculations, while the results for 
the midspreads are expected to be reliable even in the 
relatively small single-particle model space used in the 
exact calculations. 


C. Results and discussion 

In the following we present results for the level curva¬ 
ture statistics. 


1. First differential-conductance peak 

The level curvature distribution for the first 
differential-conductance peak is shown in Fig. |3] for 
four values of A/d using exact diagonalization (top 
panel) and the BCS approximation (bottom panel). The 
median and midspread are shown by the open circles in 
Fig. m as a function of A/d. 

We observe the following qualitative features in both 
exact and BCS results: (i) The level curvature distribu¬ 
tion, which is symmetric around zero at A/(5 = 0, shifts 
to negative values in the presence of pairing correlations 
and becomes asymmetric with an extended tail at nega¬ 
tive values, (ii) The modulus of the median and the dis¬ 
persion of the distribution increase monotonically with 
A/6. 

These observations are explained by the BCS approx¬ 
imation for the transition | 0 ) 7 Ve —t |0)Afe+i follows. 
The excitations that contribute to the curvature of the 
even ground state (l33l) are two-quasiparticle excitations, 
whose energies are at least 2Ae. When pairing is suffi¬ 
ciently strong, this suppresses the positive > 0) 

contribution to the observed curvature and, 

in particular, reduces the right tail of the distribution. 
For odd particle number, there are excitations that in¬ 
volve the change of the blocked orbital without changing 
the number of Cooper pairs [see Eq. dSIl) and Figs. me) 
and[2][d)]. In the limit A/6 ^ 1, the lowest excitation en¬ 
ergy for the odd grain can then be estimated as the differ¬ 
ence between two quasiparticle energies E^p+i — E^p ~ 

6'^ A'^ —A fv d^/(2A), which is much smaller than the 

noninteracting value 6. Therefore, the negative contribu¬ 
tion to the level curvature < 0) is enhanced by 



FIG. 3. The level curvature distribution P^n/do) of the first 
differential-conductance peak for A/ 5 = 0 (dotted lines), 
A/S = 0.4 (solid lines), A/5 = 1 (dashed lines), and A/5 = 3 
(dash-dotted lines). The curvature is expressed in units of the 
single-particle midspread do (the width of the middle 50% of 
the distribution). The top panel describes the results of the 
exact diagonalization method and the bottom panel describes 
the BCS results. 

pairing correlations. When mesoscopic fluctuations are 
taken into account, the lowest excitation energy fluctu¬ 
ates and can assume even smaller values, thus extending 
the left tail of the distribution. 

In the fluctuation-dominated regime A/5 < 1, pairing 
correlations are noticeably weaker in an odd grain than 
in an even grain because of the blocking effect. This can 
be observed in the exact results, which are more reliable 
than the BCS results in this regime. For A/5 = 0.4, pair¬ 
ing correlations are sufficiently strong in the even states 
|0)Ar^ so as to suppress the right tail of the distribution. 
However, in the odd states |0)Ar^+i they are not suffi¬ 
ciently strong to extend the left tail of the distribution 
or to increase the total dispersion (see the top panel of 
Fig.©. 

It is remarkable that the asymmetry effect is strong 
even in the fluctuation-dominated regime A/5 < 1, where 
the probability to observe a positive level curvature is 
small (see the exact results for A/5 = 0.4 in Fig. ©. 
We conclude that the level curvature of the first peak 
is a sensitive probe to detect pairing correlations in the 
fluctuation-dominated regime. This probe is practical 
even when only a few data points are available. 

For a single grain, level curvature statistics can be col¬ 
lected by changing the gate voltage and thus varying the 
number of electrons in the grain. An alternative way to 
observe the effect is to fix some orbital fco and study the 
curvature for the transition |0)Ar^ —>■ [4', A:oQ:o)Are-i-i as a 
function of where |4>, fcoQ;o)Afe-i-i is the lowest-energy 
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exact BCS 



FIG. 4. The median Kmed (top panels) and midspread ci„ (bot¬ 
tom panels) of the level curvature distribution versus A/5. 
The results are shown in units of do for the first (circles), sec¬ 
ond (squares), and third (diamonds) differential-conductance 
peaks, using the exact (left column) and BCS (right column) 
formalisms. The statistical errors^ are shown by vertical 
bars. The results at different values of A/5 are uncorrelated, 
i.e., they are obtained using different ensembles of random 
matrices. 


state for given ko and A/. When the Fermi level is tnned 
between a value far below ko and ko so Idf, fcoQ!o)Ar,,-i-i 
becomes the odd ground state, the total change in the 
curvature should be negative. For each value of W, the 
differential-conductance peak with the blocked orbital ko 
can be identified by its g-factor , which is independent 
of the position of the Fermi level (see Sec. IIII|) . 


2. Higher differential-conductance peaks 

We next discuss the level curvature statistics for the 
second and third differential-conductance peaks. The 
median values and midspreads are shown in Fig. |4]as a 
function of A/S by squares (second peak) and diamonds 
(third peak). 

We observe that the level curvature dispersions for the 
second and third peaks behave alike as a function of A /S 
and are larger than the level curvature dispersions for the 
first peak. However, the median values for the second 
and third peaks behave very differently from each other. 
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FIG. 5. The median level curvature ftmed (top panels) and 
midspread dn (bottom panels) in units of do as a function 
of the differential-conductance-peak number. The results are 
obtained using the BGS formalism for A/5 = 2 (left column) 
and A/5 = 5 (right column). 


To understand these results, we consider the odd excited 
states |2 )at^+i and |3)that correspond to the second 
and third peaks with energies and , respec¬ 

tively, and denote by fco the ground-state blocked orbital. 
When pairing correlations are weak, the blocked orbitals 
in the states |2)Ar^_|_i and |3)Ar^+i are usually fco -I- 1 and 
/co + 2, while the state with the blocked orbital ko — 1 is 
unresolved. However, when pairing is sufficiently strong, 
the blocked orbitals in the states |2)Ar,,-i-i and |3)Ar,,-i-i 
tend to be fco ~ 1 and fco + 1. In the limit A/5 1 and 

for an equally spaced single-particle spectrum, the ener¬ 
gies and are equal since the quasiparticle 

energies Ekg-i and Ekg+i are both V5^ -|- A^. In a re¬ 
alistic RMT-like spectrum and for finite but sufficiently 
large A/5, these energies are not equal, but tend to be 
closer to each other than to other eigenenergies. There¬ 
fore, the term in Eq. dHJ that contains the difference 
between and in the denominator is often 

the dominant term. Since such terms in the perturbation 
theory expressions have equal amplitudes but opposite 
signs for the second and third peaks, the dispersions of 
the resulting distributions behave similarly, while typical 
values of the corresponding curvatures are very different. 

Similar effect can be observed in higher differential- 
conductance peaks. In Fig. [Sj we show the medians and 
midspreads of the level curvature distributions as a func¬ 
tion of the peak number for the first nine resolved peaks 
for A/5 = 2 (left panel) and A/5 = 5 (right panel). We 
observe similar values for the level curvature dispersions 
of the 2n-th and (2n -|- l)-th peaks. However, the level 
curvature median has a negative contribution for the 2n- 
th peak and a positive contribution for the (2n -|- l)-th 
peaks. This latter effect is enhanced at larger values of 
A/5, leading to odd-even staggering in Kmed vs. conduc¬ 
tance peak number. As the peak number increases, both 
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the midspreads and medians converge to their noninter¬ 
acting (single-particle) values. 


V. CONCLUSION 

We have studied the effect of pairing correlations on 
the magnetic-field response of discrete energy levels in a 
metallic nanoparticle with spin-orbit scattering. In par¬ 
ticular, we investigated the 5 -factor and zero-field level 
curvature statistics, which parametrize, respectively, the 
first- and second-order corrections (in the magnetic field) 
to the energies measured in single-electron-tunneling 
spectroscopy experiments. 

We have shown that 5 -factors are not affected by 
pairing correlations and reduce to the 5 -factors of the 
single-particle orbitals. This conclusion follows from the 
blocking effect of pairing correlations and from consid¬ 
erations of time-reversal symmetry. It is independent 
of the strength of spin-orbit scattering and of the rel¬ 
ative effects of orbital and spin magnetisms. Thus, we 
can use 5 -factor measurements to probe the importance 
of electron-electron correlations beyond pairing interac¬ 
tions. 

In contrast, level curvatures are highly sensitive to 
pairing correlations. This can be understood qualita¬ 
tively by the dependence of the level curvature on the 
density of states, which in turn is modified by pairing 
correlations. In particular, the level curvature distribu¬ 
tion for the first differential-conductance peak, which in 
the absence of pairing correlations is symmetric around 
zero, shifts almost entirely to negative values even in 
the fluctuation-dominated regime A/d < I. In this 
regime, the pairing-induced excitation gap cannot be re¬ 
solved in a conventional spectroscopy experiment, while 
the change in the level curvature can still be observed 
and can therefore be used to probe pairing correlations. 
We have demonstrated these results by using both exact 
diagonalization and a BCS approximation. The latter 
approach can be applied in a much larger single-particle 
model space. 

The exact and BCS formalisms discussed here can be 
applied in the more general case of arbitrary spin-orbit 
scattering and in the presence of orbital magnetism. An¬ 
other possible application of the techniques developed 
here is in the calculation of the low-temperature spin sus¬ 
ceptibility of a superconducting nanoparticle with spin- 
orbit scattering. The non-divergent contributions to the 
spin susceptibility in the limit ksT <C S {ks is the Boltz¬ 
mann constant and T is temperature) resemble the per¬ 
turbative expressions for the level curvature [Eqs. (fT51) 
and (HU)]. For odd particle number, there is also a di¬ 
vergent Curie-like contribution ^ {Sl)/T, which is sup¬ 
pressed for strong spin-orbit scattering. The odd-even 
effect in the spin susceptibility of normal-metal particles 
with spin-orbit scattering was studied in Ref. l49l where 
it was shown that both the even and odd susceptibili¬ 
ties reduce to the high-temperature Pauli susceptibility 


in the limit of strong spin-orbit scattering. This corre¬ 
sponds here to the zero average curvature — Hq” 

in the noninteracting limit. The spin susceptibility of 
superconducti ng p articles with spin-orbit scattering was 
studied in Ref. but without considering the odd-even 
effect. 
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Appendix A: Many-particle matrix elements of 
in the exact formalism 

Here we derive the expressions for the many-particle 
matrix elements of between exact many-particle 
eigenstates of the even and odd grain. 

1. Even particle number 

We denote by a the time-reversed state of a (e.g., 1 = 
2 ) and by m'VJ{k} the set of orbitals containing the set m' 
and one extra orbital k that is not in m' . Using Eqs. ( 0 , 
(HI, and (HSI, and the relation 

MLm'c' = (-1)“+“'+' , (Al) 

which is valid for the phase convention ( 0 , we find 

mmn. = E 

k<.k' olol' 

X fc'y r//i~C'A Y] , kak'a') M . 

V rn'U{k'} m'U{k}J ' ’ 

m' •.k,k' ^m' 

(A2) 

This sum over Slater-determinant states is over the po¬ 
sitions of two blocked orbitals and configurations of the 
remaining Nf./2 — 1 electron pairs. Therefore, a state 
with non-zero matrix element ( 0 '|M 2 |O)Ar,, must 
have two unpaired electrons. The matrix element be¬ 
tween | 0 )Ar^ and another fully-paired state is identically 
zero because of the time-reversal symmetry. The matrix 
element of between the ground state and any state 
with two blocked orbitals 

\'^,kak'a')N^ = E] C%\rh',kak'a')N^ (AS) 

m':k,k' ^m' 

is thus given by 

{^',kak'a'\MMN. = (A4) 
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where 

s..,io.<i»i= ^ c-i. 

m' :k,k' ^m' 

, (A5) 

describes the interaction contribution to the matrix ele¬ 
ment. 


2. Odd particle number 

For a given blocked orbital fco in an odd-grain state 
m, we divide the sum in Eq. (ini into three contribu¬ 
tions 

Here consists of terms that do not contain k^, 

Mz ' consists of terms that produce the linear correc¬ 
tion (i.e., with c\^a^koa'), and consists of the re¬ 
maining terms. By analogy with the even case, the first 
contribution results in nonzero matrix elements between 
1 ^', koao)Ne+i and the states with three blocked orbitals 
(with the number of Cooper pairs reduced by one) 

(4'', kak'a'koaolMzl'^, fcoQ!o)iVe+i 

= (A7) 

Here 5^^/'k'] is defined as in Eq. (IA5I) but with C°’s 
substituted by C'^’s of Eq. dini). The phase of the matrix 
element depends on the relative positions of the three 
blocked orbitals and is not important. The second group 
of terms in Eq. (IAbI) does not contribute to the level 
curvature, while the third contribution is 


Appendix B: Details of the BCS formalism 


1. Even curvature 

The state |BCSe) defined in Eq. (l28l) is the vacuum 
with respect to the quasiparticle operators 


f i^kl — ’ 

\ Pk2 — '^k^k2 T ■ 


(Bl) 


Using the inreverse transformation, the identity 
/3fca|BCSe) = 0, and Eq. (lAlll . we find 


M^lBCSe) = ^ {ukVk' - Uk'Vk) 
k<k' 

^kl,k'lPklPk ’2 + i^kl,k'l) /3fc2^fe'l ~ ^kl,k' 2 Pkll^k'l 

+ (MfcV2)*/3L4j iBCSe), (B2) 


which leads to Eq. (ESI). 


2. Odd curvature 

Eor the odd state |BCS, fcoao) defined in Eq. (IMl) . we 
split Mz into three components according to Eq. (IA6I) . 
The effect of is similar to the effect of Mz for even 
particle number, resulting in the contribution (1551) . The 
second part - Mz - does not contribute to the level cur¬ 
vature. For the remaining part, we find [using Eqs. ((All) 
and (IMl) ] 


M^^'>\'^,koao)Ne + l = X! ^ka,koao 

koL,k^kQ 

X Y. C%\rh,ka) A^e + 1 

\7fL:k,ko^rh 

+ Y C!A\rn\{k}Ll{ko},ka)Ne+i\ ■ (AS) 

rh\k^rh,kQ^fh j 

Here the set r7i\{fc}U{A:o} is obtained from m by replacing 
k with ko- Therefore, the matrix element between two 
states with single but different blocked orbitals is given 
by 

{^",ka\Mz\-^,koao)N^+i = , 

(A9) 

where 

Dfc„fe[vI/,vi/"]= ^ c%:c^ (AlO) 

m:feo 


Mf)|BCS,A:ol) 

~ 'Y ^ka,kol'^ka(^>^ok + VkokclgicY) 

kcK,k^ko 

X n (“fcofe'+^^fcofc'4oi4o2)|vac). (B3) 

k'^k,kQ 


Therefore Mz has nonzero matrix elements between the 
states with different singly blocked orbitals. To reduce 
the computational effort, we keep only those contribu¬ 
tions that would be nonzero if the blocking effect were 
ignored, i.e., if all Ukgk’s and Vkgk’s were equal to Uk’s 
and Vk’s. For a given blocked orbital k ^ ko, there are 
two states with such matrix elements. The first is the 
lowest-energy doublet |BCS, ka) with the variational en¬ 
ergy defined in Eq. (155)) . The corresponding matrix 
element is 

(BCS,fca|M,|BCS,fcol) = 

X {ukokUkko + VkokVkko) (ukok'Ukk' + Vkgk'Vkk') ■ 

k'^k,ko 

(B4) 


and 


m if k ^ rh 

r7i\{fc} U {fco} if kGfh. 


(All) 
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The second state is the two-quasiparticle excitation on 
top of |BCS,fca) with both quasiparticles on the same 
orbital kg (so ko is not a blocked orbital) 

|BCS,fca) = cl^i-Vkko + Ukkocl^A^2) 

X {ukk' + Vkk'c\,Ak'2)\'^^'^) ■ (B5) 

k'Akyko 

The energy of such a state is given by -I- ‘i.Ekko > and 
the corresponding matrix element is 


(BCS, fca|M,|BCS, fcol) = 

X (ukokVkko — UkkoVkok) {ukok'Ukk' + Vk^k’Vkk') ■ 

k' A^iko 

(B6) 


Assuming that 

{UkQk'Ukk' H“ ^kQk' 'f^kk'^ ~ 1 ; 

k'A^Ao 


which would be an identity if the blocking effect were 
ignored, we obtain the second contribution (1571) to the 
curvature of the odd state. 


Appendix C: Differential-conductance peak heights 


Assuming a point contact between the grain and the 
electrode, the height of the differential-conductance peak 
for the transition |0)Ar^ —>■ |r2)Ar^+i in a tunneling- 
spectroscopy experiment is proportional to^ 


w'|o)^|n> = l(^^IV't(j')|0)ArJ^ + Kf2|V'j(r)|0)ArJ^ , (Cl) 

where r is the position of the contact. The field operators 
are given by a symplectic transformation [using the phase 
convention m 


^ [pk^^Aki + 9fc(i’)c^2 

k 

[-Iki^Aki +Pk{r)A 


(C2) 

(C3) 


Here [pfe(r) - q^(r)]^ and [qfc(r) Pfe(r)]^ are the single¬ 
particle spinor eigenfunctions written in a good-spin ba¬ 
sis. 


Using Eqs. and m. we find 


u;| 0 )^|O) =u;fc„(r) ^ C^C^j , (C4) 

XTrLiko^rh j 

where 

WfeoW = bfeo(i-)P + |g/co(i-)P ■ (C5) 

The peak height is now conveniently expressed as the 
product of a single-particle and many-body contribu¬ 
tions. The single-particle contribution Wk^ (r) is simply 
the probability density for an electron in the orbital /cq- 
The many-body contribution is generally nonzero for any 
eigenstate with a single blocked orbital. In the noninter¬ 
acting limit, it is nonzero only for states |H) 7 v^+i shown 
in Fig. Iljc), for which it is equal to 1 and the total peak 
height reduces to Wk^ (r). 

In the BCS formalism, the even ground state is given 
by Eq. (l28l) . and any odd state with a single blocked 
orbital fco (i.e., not only the lowest-energy state) can be 
written as 

\^)n,+i = (uk + ikcAL) |vac). (C 6 ) 

k^ko 

For these states, we obtain 

u;|o)^|n) = Wko (r) ul^ {ukUk + VkVkf ■ (C7) 

k^ko 

In the limit A/6 ^ 1, the blocking effect is negligible 
and all the odd states are the quasiparticle excitations 
built on top of |BCSe). For the one-quasiparticle excita¬ 
tion |BCSe), Uk = Uk and Vk = Vk, and the peak 
height re(3uces to 

w\o)^\n) =Wko{r)ul^. (C 8 ) 

Any other odd state (IC 6 I) with the same fco can be writ¬ 
ten as PAoAAk' 2 ---\^^^^) fo’' ^ ^ 0 - The 

product P\ii 0 ki 2 results in Uk' = —Vk' and Vk' = Uk', 
which yields ui|o)^|a> = 0 . 

We conclude that, in the limit A/S ^ 1, a final state 
must have the lowest energy among the states with the 
same blocked orbital. The interaction contribution to 
the peak height (ICSp is given by ; therefore, the peak 
height for such a state is suppressed when the blocked 
orbital is far below the Fermi level and reduces to the 
noninteracting value when it is far above the Fermi level. 

In the RMT framework, the single-particle contribu¬ 
tion (IC5I) is derived from a component of the spinor eigen¬ 
vector representing the state |fcoao)i. It is not important 
which component is chosen. In any of the Gaussian en¬ 
sembles, the probability for Wka (r) to be small is expo¬ 
nentially suppressed. 
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